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Abstract 
as 

Q^l We study the correlation between e'/e and A/ = 1/2 rule in the frame- 

o ' work of non-linear a model including scalar mesons. Using this model we 

^ ■ estimate the chiral corrections by resumming to all orders in Chiral Perturba- 

tion Theory, the contribution of a class of diagrams within the factorization 
approximation. With these matrix elements and changing the scalar meson 
kS ' mass, we find that there is correlation between e'/e and AI = 1/2 ampli- 

C^ • tude. However, it is difficult to explain both e'/e and A/ = 1/2 amplitude 

simultaneously. In order to be compatible with e'/e, typically, about half of 
A/ = 1/2 amplitude can be explained at most. Our result suggests there 
may be substantial non-factorizable contribution to CP conserving K -^ tttt 
amplitudes. 
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I. INTRODUCTION 

According to recent measurements of direct CP violation in K ^ tttt decays, e'/e is 
O(10~'^) |jl|,@]. Theoretical prediction ranges around 10~^ ~ 10~^. It strongly depends on 
the hadronic matrix elements of QCD P| and EW penguin operators @-§|. 

An interesting possibility is suggested as an explanation for large e'/e = O{10~^) in the 
standard model 0. The authors argue that the mechanism which enhances A/ = 1/2 ampli- 
tude may also enhance e'/e and it naturally leads to the measured values. The enhancement 
comes from the Feynman diagram which includes the scalar a meson as intermediate state. 
The mechanism was found in the framework of linear a model |^,^. Though qualitative 
picture of the linear a model may be correct, for quantitative analysis, some improvement 
can be made. In the linear a model employed in Ref. [^, a meson mass can be written in 
terms of the physical quantities Fk, -Ftt, Mk and M^ and its mass is predicted to be about 
900 MeV. SU(3) breaking ratio irtu/ms is also determined by the same input and numer- 
ically it is around 1/30. The enhancement factor of a QCD penguin operator is written 
as FkH^F^ - 2Fk) = {Mna^ - M^^)/{M^^ - Mk'^) ■ (Fk^/F^^) ~ 2. These relations and 
numbers are specific predictions of the linear a model. Because the dynamical property of 
the linear a model is not same as that of QCD, these relations and numbers may be taken 
as semi-quantitative ^. 

In this paper, we study the matrix elements of the QCD and EW penguins with the non- 
linear a model including scalar mesons. The model is built with chiral symmetry as a guide. 
It is more general than the linear a model and less dependent on dynamical assumption. 
The cost is that it has more parameters. They can be determined with the experimental 
measured quantities, i.e., decay width, mass spectrum, etc. Still a meson mass is left as a 
free parameter because the spectroscopy of a meson /o(400 — 1200) allows wide range for 
the mass. (See Refs. ||T0Hl5[ for scalar meson mass spectroscopy.) We study how QCD and 



EW penguin matrix elements depend on the mass of a meson. 

Here we write a few words on the difference between our approach and the conventional 
treatment of the resonances in Chiral Perturbation Theory (CHPT). CHPT is a systematic 
treatment, in the sense of the small momentum expansion, and describes low energy processes 
involved pions and kaons well. However higher order terms in CHPT are of great importance 
if threshold of the a meson is near to kaon mass. Thus it is very interesting to explore non- 
linear a model with scalar resonances in kaon decays. Correspondence between non-linear 
a model with resonances and CHPT was well discussed in Refs. [^|l7l at order of p^ and 



their results indicate that the resonance contributions dominate the low energy coupling 
constants in the strong part of p^ chiral lagrangian. In our approach, we compute the 
full contribution of the a meson to the factorizable part of QCD and EW penguins using 
non-linear a model with scalar resonances, so that a certain class of higher order terms of 
p"'{n > 4) in momentum expansion are included. We show how these higher order terms 
contribute to observable in kaon decays. 

As a phenomenological application, we compute e'/e and A/ = 1/2 amplitude in isospin 
limit. We use the Wilson coefficients for four-fermi interaction at fi = 0.8 ~ 1.2 GeV in the 
next to leading log (NLL) approximation [|18H20[|. The four-fermi operator is factorized into 
products of color singlet currents (or densities) and they are identified with those of the a 
model [^. The density x density type operators are enhanced for small strange quark mass 



by a factor of (l/rris)^- Therefore numerical values for the strange quark mass are important. 
As for the strange quark mass, we choose the range which is suggested by QCD sum rule 
p^ and lattice simulations |2^. We also study the correlation between e'/e and A/ = 1/2 



amplitude. This is done by varying the a mass, strange quark mass nig, and factorization 
scale /i. By studying the dependence of e'/e and A/ =1/2 amplitude on the a meson mass, 
we search for the range of the mass which may reproduce both e'/e and A/ = 1/2 amplitude. 
The paper organized as follows: In section II, we summarize the outline of the computa- 
tion e'/e and AJ = 1/2 amplitude. In section III, we derive the matrix elements of penguin 
operators. In section IV, numerical results of e'/e and A/ = 1/2 are summarized. In section 
V, we discuss the implication of our results. Some useful formulae are collected in appendix. 



II. A/ = 1/2 RULE AND e'/e IN THE STANDARD MODEL 



In this section, we summarize our notations and show outline of computation of e'/e 
and AJ = 1/2 amplitude. Some details of definitions of isospin amplitudes can be found in 
appendix A. We start with the effective hamiltonian for AS" = 1 non-leptonic decays |T^, 
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Heff = -^VudV*, J2{zi + ryijQi + h.c, 
v2 i=i 



(1) 



where r = —(ytdV^*)/(yudV*g). The isospin amplitudes oi K ^ rni are defined as 
{I\Heff\K^) = iaiexpi6i, {I\Heff\K^) = —iai*expiSi. e' is expressed in terms of qq 
and 02. 

f 1 Rea2 ,Ima2 Jmao, 



exp{i{52 - So 



\pl Reao Rea2 Reao 

In the factorization approximation, aj s are written as, 

Gf 
Reai = -j=ReVudV*, 



-)}• 
2^^ 
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[z2^-l+"f){I\Q2^-l\K') + iz2^+'-j^ 



){I\Q2^\K') 



'V2 



X 
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ImVtdV,l 



y2 



y2i~i^ 



E (?/2-l + W)(^\Q2^-1\K') + fe + ^){I\Q2^\K') 



N, 



Nr 



(3) 



(4) 



where the matrix elements {I\Qi\K^) is defined in large Nc limit. As we discuss in detail 
in the next section, we compute the hadronic matrix element in large Nc limit, i.e., we 
factorize the four-fermi operators into products of color singlet currents (densities). The 
currents (densities) are identified with those of the chiral lagrangian. The factorization 
scale is chosen at 0.8 — 1.2 GeV, i.e., below charm quark mass rric. In the factorization 
approximation, this choice is mandatory because above rric, real part of the Wilson coefficient 
of QCD penguin operators is zero and it is born below rric due to the incomplete cancellation 



of GIM mechanism. About Wilson coefficients, we use NLL approximation |T8|-pO| and 
compute them at the factorization scale. Combining the matrix elements with the Wilson 
coefficients, a/ s are given as. 



Reao -- 
Rea2 -- 
ImaQ 



GfX 



2z2 + 3^4 - zi + 3ze 



V3 



\—[zi + Z2] 



'^^^^-{AX'fXv 



V2 



Ima2 = r(^A^) Xt] 



X 



1/4 + 



y7-y9 + yio 



Yf 



6 , 2 Ys — -^Yq 



^6^+3^/8 ^ 



2 Ys 

-1/7 + 1/9 + 1/10+ 3^/8^ 



(5) 



where matrix elements are denoted by X, Yq, and Yg. X is the matrix element of 
current X current type operators, Ig corresponds to the matrix element of a density x density 
QCD penguin operator, Kg is the matrix element of EW penguin operator. Their derivation 
and precise definition will be given in the next section and are summarized in Table I and 
II. 



III. NON-LINEAR a MODEL INCLUDING SCALAR MESONS AND THE 
MATRIX ELEMENTS OF QCD AND EW PENGUIN OPERATORS 

The non-linear a model with higher resonances are studied in [T^|T^ . In K ^ tttt decays, 
in large N^ limit, scalar meson may contribute to the matrix elements of density x density 
type four-fermi operators {Qq, Qs)- For current x current type four-fermi interactions, the 
amplitude is proportional to the form factor of semi-leptonic decay, i.e., /^{Mk'^ — M^^) + 
f-M^. Because the form factors f±{q^) near soft-pion limit (g^ = M^^) are important, vector 
mesons contribution to the form factors is small and their effect can be safely neglected. 
Therefore we include only scalar mesons in chiral lagrangian. 



f 



91, 



C = ^—TrdUdU^ + B TrM{U + U^) + —TrdUdU^iSi'^ 
+ g2 TrM{iSi + S}SS}) + Tr{DSDS - KCS^), 

where A^ = diag{mu, riid, rris), U = exp {i27i/ f) = C,"^ and 5 is a scalar nonet field, 

/a + 6^ V26^ V2k+ 

S = -\ x/25° a -6^ V2^' 

\ V2k- V2k,o V26ss 

DS is covariant derivative and is defined by. 



(6) 



(7) 



D,S 



«//p 



d^S + i[a//^,S], 

2i 2iP 



+ ... 



(8) 
(9) 



In the lagrangian, the scalar mesons couple to pions through two terms denoted by gi and §2- 
One is a coupling in SU(3) limit and the other is a coupling with SU(3) breaking. The mass 
splitting term for the scalar nonets and isospin breaking effect are neglected. By shifting 
the scalar meson fields from their vacuum expectation value, 



S^S+<S>, <S>-- 
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M, 



(10) 



we obtain the mass formulae and decay constants [T^. They are given in the appendix 
B. The parameters in the chiral lagrangian can be written in terms of physical quantity 

Fk, -Ftt, Mk, Mt^ and quark masses mu{= irid) and mg. 
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+ 



2(1 + A) 



F1M| 



4A 2(1 + A) 



;i-A) 



m. 



(11) 
(12) 
(13) 



where A = rriu/ms. For computation of the weak matrix elements, we need strong interaction 
vertices. They can be found in appendix C 

In our calculation using the lagrangian Eq.(^, a certain class of the higher order terms 
in the CHPT are summed up due to the effect of the scalar resonance exchange. These are 
very important in the process K ^ tttt if the a meson mass is light as kaon mass, M„ ~ M^- 
Though systematic treatment of momentum expansion in the CHPT is lost, a class of the 
higher order terms in the CHPT are automatically summed up using lagrangian of Eq. (^ . 

Now we turn to the matrix element of QCD and EW penguins. The explicit derivation 
is given for two density x density type operators Qq and Qs- Their definitions are. 



Qe = -8 Yl (siqRliqRdL), 

q=u,d,s 

Qs = -12 Yl {sLQRieqiqRdL) 

q=u,d,s 
Qs = -S{sLUR){uRdL), 



Q, 



2^6, 



(14) 



(15) 
(16) 



where the subsidiary operator Qs is introduced. These operators can be written in terms of 
meson fields by identifying quark bilinear as corresponding density. 



QRQi 



-BUi. 



After some algebra, we express Qq in terms of the meson fields. 



(17) 
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i glirus - ruu) L^ ^|(m, + m„) \ K_ 
MJ [ MJ / Fi 
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g^[ms -rriu] 
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M„ 
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K^TT^^ 



M. 



2B + gi 



2 [ms + niu) 



71^ K 



M/ 



+ 25 + 



FkF^ 
gl2mA aK^' 



M/ 



'K 



(18) 



There are four diagrams which may contribute to K'^ -^ tt^tt*' amphtude Yq. (See Feynman 
diagrams in Fig. 3 - Fig.5) They are classified as follows. 

1) The diagram in which K^ decays into K'^tt^tt^ through the strong vertex and subsequently 
K^ vanishes into vacuum through the Qq ( Tfadpoie)- 

2) The diagram in which K^ directly decays into 27i^ {Tdirect}- 

3) The diagram in which K^ is converted into a and subsequently a decays into 2-n^ (Tfj-poie)- 

4) The diagram in which K decays into k, vr*^ through strong vertex and k is converted into 

'' y-'-K-pole)- 

The sum of the contribution is denoted by Yq and it can be simplified as, 
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2FIM,' 
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(19) 



where R = Fk/ Ft, and 5k{tt) = Mk{-k)IM^. We also introduce the auxiliary quantities, 



B = B 



2M„ 



Bo = B + 






(20) 



They can be written in terms of physical quantities. 



B 



2m,(l + A)^ 



5n 



MlFl 
4m„A' 



(21) 



The matrix element of EW penguin operator Qs is straightforward. Technically we split Qq 
from Qs so that we do not have to repeat the calculation of Qq. The rest is called Qs and 
given by. 



Qg 



■. 71-^+ .g2Bo{Bo + B)K'7r+iT- 



V2 



+ ■ 



(22) 



In K^ —>■ TT+TT amplitude, there are two contributions to the hadron matrix element of Qs, 
i.e., (a) K-pole contribution and (b) direct contribution. The sum is called Fg and is given 
by. 



Y. 



Bo 



V2FkF^ M^^ - Ml [ 
+26v^-^5o (b + 5, 



gig2Ml - 2gl \ (m, + 3m) 



m„ — m) 



Mj 



(23) 



Keeping leading terms of 1/M„ the matrix elements 1^6; ^8 reduce to the well known results 
Yq,Y^ [|^, which correspond to those in the leading order of momentum expansion and in 
large Nc limit. 



2 \ 2 / /\#2 



Y,' = -AV2z {Fk - F^) ( ^^ \ , Y^ = ^^iF, 



M 



K 



TTls + mu, 



(24) 



This approximation is valid only when M^ -C Mq-. In the next section we will show how 
the values of the matrix elements of the density x density operators are different from Yg°, Y^ 
numerically. 

The matrix elements of the other current x current operators Qi{i ^ 6, 8) are also shown 
in Table I and II, and are expressed by a single amplitude X: 

X = iV2f{MK^-M^^). (25) 



IV. NUMERICAL RESULTS 

In this section, we first estimate the hadronic parameters -06,8 corresponding to the 
matrix element Qq and Qs in the factorization approximation. We compare our results with 
those from the linear a model. As an application, we also compute e'/e and ReaQ, Rea2. 
This is done in isospin limit. 

The conventional bag-factors Bq , B^ , which people often refer, is defined by the 



following equation in our notation, 

r1/2 _ X^ o3/2 _ 2^8-^ .2fi^ 

As explained in section III, using M,r, Mk, F.„, Fk as inputs our model can be described 
by three free parameters in the lagrangian M^, iris, A = niu/ms and another parameter ix 
in the matching process which is the factorization scale. 

We find that our model predicts that the factorizable part of -Bg ranges around 1.6 ~ 3.0 
depending on the a meson mass M^.. The quark mass dependence is not significant. On the 
other hand, B^ ranges around 0.7 ~ 1.1 depending on A. We varied A in the range of 
1/20 — 1/30 and smaller A gives larger value of B^ . M^ dependence is negligible for B^ . 
The numerical value is given in Table III. 

Let us now compare our results with those from linear a model 0. The factorizable 
bag-factors are given by the following equations. 

Bl^' = ^ 2, (27) 

^ 3F^-2Fk ^ ' 

bT ^ §^ ~ 1.2, (28) 

Because there are only four parameters in the linear a model lagrangian, after using 
Mt^^x,Ft^,k there are no free parameters left, so that the model predicts A = 1/30 and 
M<,'~0.9GeV. 

3/2 

From Table III we find that Bg from our model with A = 1/30 and that from the linear 

1/2 

a model are consistent. It can be seen that the factorizable part of Bq for Mo- = 0.9 GeV 
is around 1.5, which is smaller than the linear a model result. 



Next we apply our result to e'/e and Reao, Rea2. For numerical computation of e'/e, we 
use the experimental values for Reaj. 

We have calculated the next to leading order Wilson coefficients in NDR scheme, we 



could reproduce the numerical values tabulated in Ref. [|^ to a good extent. We chose the 
following values for the computation, nit = 165.00 GeV, niw = 80.20 GeV, nib = 4.40 GeV, 
rUc = 1.30 GeV, l/aqED = 129.0, sin'^{9w) = 0.230, AgJ.^ = 0.226 GeV, AgJ.^ = 0.325 
GeV, a^^(m^)(5) = 0.11799. We list the Wilson coefficients Zi,yi at scales /i = 1.2, 1.0, 
0.8 GeV in Table IV. In this calculation, we used the anomalous dimensions at the NLO by 
Buras et al. 

As was explained in section II, in the leading order in large Nc expansion, Reaj, Imaj 
are obtained by multiplying Wilson coefficients Zi{^),yi{fi) with the matrix elements of 
Qi{fj,), ■ ■ ■,Qio(/i) in our model, where /i is the factorization scale which is assumed to be 
0.8 - 1.2 GeV. 

Here we should remark one point about the quark mass. In large N^ limit, we approximate 
the matrix elements with Qq{^) operator by the product of matrix elements with scalar 
quark operator at scale /i. Using PCAC relation, we then convert them to F^M'^/'ms^n). 
Here, the scale of the strange quark mass should be the same scale /i. Therefore, when we 
substitute the mass parameter m^ in our final result, we should run the quark mass to the 
factorization scale /i=0.8,1.0,1.2 GeV. For example, ms(2GeV)= 80 — 120 MeV corresponds 
to m,(0.8GeV) = 136 - 204 MeV. 

Figure ^ shows the dependence of t] from our model on the scalar resonance mass 
M„. Here, we take mf^(2GeV)=80,120,180 MeV, which cover the recent QCD calcula- 



tions P2|,p3, and the scale jjl is chosen to be 0.8 GeV. Upper and lower lines correspond to 



the maximum and minimum values of e'/e, respectively. We find that when Mo- is larger 
than 0.8 GeV in order for r] to lie within 0.27-0.52, which is favored by other measurements 
of CKM parameters, ?7i*^'^(2GeV) should take rather small value 0.09 — 0.12 GeV. These 



values are consistent with recent lattice QCD calculations |2^ but smaller compared with 



QCD sum rule results ||22|. On the other hand, as M^ becomes smaller the Qq amplitude 
gets enhanced, in order for r] to lie within 0.27-0.52, larger value of ?Ti^'^(2GeV) is preferred. 

Finally, in Figure 0, we show the correlation of ao/oQ^'' and e'/e by changing the scalar 
meson mass from 0.6 GeV to 1 GeV. We take three different values of m^^(2GeV), which 
are 80, 100 and 120 MeV. We take three values for the factorization scale /x, which are 0.8, 
1.0 and 1.2 GeV. We used the experimental values of Reai for the numerical analysis of 
e'/e. CP violation parameter is chosen to be 0.3 in the figure. The shaded region is the 
experimental data from KTeV and NA48 for e'/e at 2-cr confidence level. 

We find that e'/e can be easily explained in our model by a suitable choice of the pa- 
rameters. Typical value of Mo-, and m^^(2GeV) is around 0.8 GeV and around 80 MeV 
respectively almost independent of the factorization scale /i. On the other hand, Rea^ in 
our model are smaller than experiment. We find that it is quite sensitive to the factorization 
scale /i, and as /i gets smaller, Reao becomes larger towards the experimental value. For 
/i=0.8 GeV, aQJa^Q^^ is around 0.5 — 0.6. 

The sensitivity of Oq amplitude on n can be understood as follows. The Wilson coefficient 
zq{^) vanishes when the GIM cancellation between the charm penguin and up penguin loop 
is exact. In our calculation, since we take the MS scheme, the cancellation is exact above 
charm threshold. Therefore, Zq{^) takes nonzero value only when /z < iric. Since the 
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factorization scale is very close to the charm threshold, the result of zq{^) changes quite a 
lot. 

In contrast, the Wilson coefficient yeifi) vanishes when the GIM cancellation between 
the top penguin charm penguin loop is exact. Since top decouples already below M^, 
this cancellation is completely violated and yeln) takes nonzero value from the start and 
keeps growing all the way down to the factorization scale. Since /o(7(Mvy/1.2GeV) and 
log{Mw/0.8GeV) is almost identical, ye^fi) is not so sensitive to the factorization scale. 
About Rea2, our result is about 1.5 times larger than the experimental value. 

V. SUMMARY AND DISCUSSION 

In this paper, we study the correlation of A/ = 1/2 amplitude and e'/e in the framework 
of non-linear a model including the scalar mesons. We have calculated the matrix elements 
of the QCD and EW penguins using that model and within the factorization approximation. 

We can not find the scalar meson mass region which is compatible with both e'/e and 
A/ = 1/2 amplitude simultaneously. The reason is follows. We can read from Fig.2, the 
maximum allowed value for {e' /e)/r] is about 0.01. The bag-factor -Be^ required for e'/e 
is at most 2 — 3, which corresponds to M^\min — 0.6 — 0.7 GeV. In the range of the scalar 
meson mass, about half of A/ = 1/2 amplitude may be explained. Therefore, if we impose 
the e'/e constraint, we can not explain the whole A/ = 1/2 amplitude. Moreover, e'/e is 
rather stable for the change of the factorization scale. This suggests that the prediction of 
e'/e may be more reliable. Though there is strong correlation between AJ = 1/2 amplitude 
and e'/e, we conclude the understanding of AJ = 1/2 rule may not be complete. 

Finally, we argue what kind of effects may remedy the problem. Because QCD penguin 
Qq is born just below rric, the coefficient is not stable about the change of the factorization 
scale around 1 GeV. In the scheme, in which GIM cancellation is incomplete above /i > rric, 
the leading order results of the Wilson coefficient of Qe becomes larger by a factor of 2 [^ . 



This effect was not incorporated in the Wilson coefficients of NLL approximation employed 
here. Therefore the same effect may further enhance the Wilson coefficient of Qq used in our 
analysis. We also note that the real part of the A J = 3/2 amplitude is larger by a factor of 1.5 
than the experimental value. This may tell us that there is some suppression (enhancement) 
coming from the low energy evolution (pion loops) for CP conserving A/ = 3/2 (A/ = 1/2) 
amplitudes [p^-[29[|. A plausible explanation is that the non-factorizable contributions are 
very large. Including these effects may help for the entire understanding of both e'/e and 
A/ = 1/2 rule. 
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APPENDIX A: AI = 1/2 RULE AND e'/e 

Here we summarize isospin amplitudes and e'/e. 

T[Ks -^ 7r°7r°] = Br{K, -^ tt^tt") x l/r„ (Al) 

T[Ks -^ TT+TT-] = Br{Ks -^ TT+TT-) X l/r„ (A2) 

r[i^+ ^ 7r°7r+] = Br{K+ -^ 7r°7r+) x l/r+, (A3) 

where 

r, = (0.8927 ± 0.0009)10~^°sec, (A4) 

r+ = (1.2386 ± 0.0024) 10-^sec, (A5) 

— = 5.3142 X lO^^^MeV, (A6) 

Br{Ks -^ 7r°7r°) = 31.39 ± 0.28%, (A7) 

Br{Ks -^ TT+TT") = 68.61 ± 0.28%, (A8) 

5r(fs:+ ^7r°7r+) =21.16 ±0.14%, (A9) 

|/ = 0) = ^|7rV) + ^|7r+7r-), (AlO) 

|/ = 2, J3 = 0) = -yf IvrV) + ^k+TT-), (All) 

|J = 2,/3 = l) = |7r°7r+), (A12) 
where iTf^vr^) and |7r*^7r+) are the symmetrized states defined as 

\7l + 7l-) = (|7r + ) X \7l-) + \7l-) X \7l + ))/V2, (Al3) 

|7r+7r°) = (|7r+) x |7r°) + |7r°) x \n+))/V2. (A14) 
We can write the decay rates in terms of the isospin amphtudes: 

{7r^n'^\Hu,\Ks) = iV2{\ -Reaoexpi5o - \ -Rea2expi52}, (A15) 

{-K^-K^\Hyj\Ks) = iV2{\ -ReaQexpi6o + \ -Rea2eyipi62}, (A16) 

{7r+7r'\H^\K+)=t^h2, (A17) 

where ia/exp(i(5/) = {I\Hu,\K°), I = 0,2 and \Ks) ~ ^jd^o) - l^o))- With the definition, 
we can write: 

2 1 

T{Ks -^ 7r°7r°) = P-\Reaoexpi6o - V2i?ea2 exp i^s ^ - , (A18) 

o ^ 

2 1 

r{Ks -^ tt+tt") = P-\V2Reaoexpi6o + Rea2expi62\'^-, (A19) 

T{K+ -^ 7r°7r+) = P^i?eao', (A20) 

10 



where imaginary parts are neglected. P is a phase space factor of two body decay and is 
defined as: 

P = i Jl - AMl/Mk^ (A21) 

= 3.34919 X 10-^{MeV'). (A22) 

Here we use M^+ = 493.677 MeV and M^ = (M^o + M^+)/2 = 137.273 MeV. With these 
definitions, we obtain, 



Br{K+ -^ 7r07r+) r, 3 Rea-, 



2 



Br{Ks -^ ttOttO) + Br{Ks -^ tt+tt-) r+ 4 Rea2^ + Reao 
We can extract the following ratio and values for ao and 02, 

Reao 



2- 



(A23) 



22.15, (A24) 
Rea2 

Rea2 = 2.114 x lO-^(MeV), (A25) 

Reao = 4.686 x lO-^(MeV). (A26) 



APPENDIX B: DECAY CONSTANTS, MASS FORMULAE 

In this appendix, we collect the formulae for the decay constants and masses which can 
be derived using Eq.(6). 



1 r. R , .92^u'^ 

—-\4:Bmu + 4: o h 



M^ = —[ABm^ + A'-i-^], (Bl) 



r2 1 roD/ , \ , gKjns + rriuf 



- ^:[1 + ?^], (B3) 



^^ 2Ma''f 



Fk = ^=[1 + ''''Z:2'2 ]^ (B4) 



^. = 1 + ^^, (B5) 



2M.V' 



^^ = 1 + o.^ 2.2 ' *^B6) 



where Zyr and Zk are wave function renormalization constants, F^ is 92.42 MeV. 

APPENDIX C: LAGRANGIAN 

Here we record the part of the lagrangian which is relevant for calculation. 
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£ 



vr 



^dK'^dK'^ irOirOSTrO" n'^dn^K'^dK'^ + K'^dK'^) f 1 1 



4-71 



12F| 






24 



i^l ^ F^^ 



s '^u 



(CI) 



TT-' K 



r - ^ 



oai^o^ 



0" 



-^2 



a 



/(7r0)2 



i^l 



k"^^^: h /i.e. 



F.F, 



■E^K 



rriu + 



K^K" frriu + nis 



p2 ■■"« 

IT 



P2 



K 






+ /i.c. 



(C2) 



where SCs-k-k and ^£47^ come from the covariant derivative term. (See Eq.(8) and Eq.(9).) 
Their exphcit forms are, 



SC 



g2{ms - mu) 



'^^ M^HFkF^ 



SC 



4tt 



glirus - m^y 



dR°{KdTT°) + h.c. 



(C3) 
(C4) 



APPENDIX D: THE MATRIX ELEMENT OF Qa 
We give the derivation of the matrix element of Qq. 

-'6 -'■tadpole ~r J- direct ~r J- a-pole ~r J- a- tad ~r J- n-pole ~r J-K-tad- 

The exphcit expression of the parts of Eg . ([191) is given by, 



-^ tadpole ~r -t direct \ A% 



M/Fk 



[nis - mjB 



(Dl) 



\fI^ F^ M^FIF^ FlF^M^^Ml^ F^F^M^'MI 



^ a-pole ~r J- a-tad 
-^ K-pole ~r -L K-tad 



iV2B 



Fkf:^{mi - M,'') 

iV2B 

FkF^{M,''-MI) 



{-gig2 {2Ml - Ml) + Sglm^] 



_ g'^jrus - mu){ms + rriu) glinis - niuf 



MlFlMl 



2FIMIMJ' 



Ml 



2MI ,,_ ^,o(. Mr 



X { -gig2Mj^ + Ag^ims + m„) - 2(m, - m„)^2TyT " 2("^s - "^«)^2 U ^ TtI 



M. 



M2 



(D2) 
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FIG. 1. Allowed regions for M^-rj from e'/e data at 2-cr confidence level. They are shown for 
three different values for the strange quark mass, i.e.,?7is(2GeV)= 80, 120, 180 MeV. 
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FIG. 2. Correlation of ao and e'/e. The stars on the lines correspond to Mq- =0.6, 0.7, 0.8, 0.9, 
1.0 GeV, ms(2GeV)=80, 100, 120 MeV at factrization scale /i =0.8, 1.0, 1.2 GeV. 



16 



71 



71 



K. 



> ® 



71 




FIG. 3. Feynman diagrams for Tdirect + Ttadpole- 



n 



71 




FIG. 4. Feynman diagrams for T„-poie + T^-tad- 
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FIG. 5. Feynman diagrams for T^-pole + T^-tad- 



17 



TABLES 



vr^vr^ |7r+ > x|7r > |7r°>x|7r+> 

"Q^ ^X liR^ 

Q2 X X/V2 

Qs 

Qa X X 

Q5 

Qe Yq Ye 

Qj 3X/2 -'iX/2^/2 

Qs yg ^8/^2 



Qg -3X/2 3X/2\/2 
Qio -X/2 X X/V2 

On X X/V2 

TABLE L The matrix elements, where X = i\/2f{Mk — M^^ ). Ig and Yq are defined in the 
text. 



Go 02 



Q2 75X ^X 

Q3 

Q4 \/3X 

Q5 

Qe VSYe 

Qs Vs - ^Y, ./fyg 



73^8-2^6 V3-^8 

Qg 2~^ V 2-^ 



2 ^'- V 2 

TABLE n. Contribution to isospin amplitudes. 



A = 1/20 


M, (GeV) 


0.55 


0.6 


0.7 


0.8 


0.9 


1.0 


-De 

„3/2 
-^8 


5.27 
0.70 


3.29 
0.73 


2.22 
0.77 


1.84 
0.79 


1.64 
0.81 


1.52 
0.82 


A = 1/25 


M^ (GeV) 


0.55 


0.6 


0.7 


0.8 


0.9 


1.0 


-De 

„3/2 

-DS 


4.81 
0.90 


3.06 
0.93 


2.12 
0.96 


1.78 
0.99 


1.61 
1.00 


1.50 
1.01 


A = 1/30 


M, (GeV) 


0.55 


0.6 


0.7 


0.8 


0.9 


1.0 


-fc'e 

„3/2 
-^8 


4.48 
1.13 


2.90 
1.15 


2.05 
1.17 


1.75 
1.19 


1.60 
1.20 


1.51 
1.21 



TABLE III. Bag-factor 



Wilson coeff. 


fi = 80.2 GeV 


fi = 1.2 GeV 


H = 1.0 GeV 


11 = 0.8 GeV 


yi 


0.0 


0.0 


0.0 


0.0 


2/2 


0.0 


0.0 


0.0 


0.0 


ys 


0.0014715 


0.03058 


0.03335 


0.03722 


2/4 


-0.0019375 


-0.05871 


-0.05884 


-0.05844 


2/5 


0.0006458 


0.00311 


-0.00168 


-0.01384 


2/6 


-0.0019375 


-0.09797 


-0.11672 


-0.16226 


yj/aQED 


0.1262367 


-0.03714 


-0.03822 


-0.04038 


y&/0LQED 


0.0 


0.14352 


0.17174 


0.23136 


y<i/oiQED 


-1.0606455 


-1.46549 


-1.54058 


-1.69377 


yw/oiQED 


0.9 


0.57829 


0.68795 


0.89882 


Zl 


0.0526643 


-0.45108 


-0.52381 


-0.64505 


Z2 


0.9812457 


1.23913 


1.28816 


1.37464 


z-i 


0.0 


0.00674 


0.01353 


0.03059 


za 


0.0 


-0.01980 


-0.03704 


-0.07439 


Zb 


0.0 


0.00569 


0.00784 


0.00844 


z% 


0.0 


-0.01950 


-0.03698 


-0.08023 


Zj/OLQED 


0.0 


0.00940 


0.01249 


0.01989 


Zs/OLQED 


0.0 


0.00349 


0.01551 


0.04725 


Zg/OLQED 


0.0 


0.01127 


0.02019 


0.03993 


Zia/aQED 


0.0 


-0.00219 


-0.00893 


-0.02287 



TABLE IV. List of Wilson coefficients yi,Zi. 
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